Scalability of GHZ and random-state entanglement in the presence of decoherence 
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We derive analytical upper bounds for the entanglement of generalized Greenberger-Horne-Zeilinger states 
coupled to locally depolarizing and dephasing environments, and for local thermal baths of arbitrary tempera- 
ture. These bounds apply for any convex quantifier of entanglement, and exponential entanglement decay with 
the number of constituent particles is found. The bounds are tight for depolarizing and dephasing channels. We 
also show that randomly generated initial states tend to violate these bounds, and that this discrepancy grows 
with the number of particles. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Yz 



I. INTRODUCTION 

Impressive experimental progress in the manipulation of 
composite quantum systems towards real-world applications 
of quantum information and communication theory has taken 
place over the last few years. The controlled production of 
genuinely multipartite entangled states plays a central role in 
this program, and is crucial for the scaling of existing toy 
quantum protocols to mature technologies. Among those, the 
W HI and Greenberger-Horne-Zeilinger (GHZ) [2| types of 
entanglement play a paradigmatic role, since they incarnate 
characteristic traits and subtleties of multipartite entanglement 
which allow, e.g., to implement protocols for secure quantum 
communication. 

Three photon W-type entanglement HI, and five- and 
six- photon GHZ entangled states [5, 6| have already been 
observed. Very recently, a ten party GHZ state was pro- 
duced using five hyperentangled photons 0. GHZ states 
have also been reported in cavity QED experiments (8], and 
for three |9), four ifTOl and six ifTTI trapped ions. Eight ions 
were prepared in a W state 1121 . and the controlled genera- 
tion of different multipartite entanglement families was shown 
in fl3l . Furthermore, the recent implementation of a robust 
and extremely-high-fidelity entangling gate of two ions lfl4l 
opens the way to the controlled production of GHZ states of a 
few tens of ions. 

However, it is known that scaling multipartite entangled 
states up to many constituents is haunted by the decoherence 
processes arising from the unavoidable and (most times) detri- 
mental interaction of the system degrees of freedom with the 
environment. Complete disentanglement may even occur at 
finite times 03J [16] [17] HU [H Ej] ED [22] [26] [27] IH, as ex- 
perimentally demonstrated in [23, 24] [25). It is therefore vital 
for experimental implementations to predict the time scales on 
which appreciable amounts of entanglement will prevail, un- 
der realistic assumptions on the environment coupling. Theo- 
retical studies of disentanglement dynamics for large systems 
were presented in Q6] [Sl[26] W\ [28] [29] [30] HQ El [33]|3l, in 



particular for GHZ and W type states. While Refs. 
considered the scaling of the disentanglement time with the 
system size, Refs. 1261 [3T1 [32l put the focus on the scaling 
of the short time scale behaviour of entanglement, what al- 
lows for a much better characterization of the robustness of 
the initial state's entanglement under noise. In particular, it 
was shown in PD that even if the disentanglement time of 
GHZ states grows with the number of parties N, the residual 
entanglement (there quantified by the negativity [35 1) is re- 
duced to arbitrarily small values at times which decrease with 
N. 



In the present paper, we expand the above studies to larger 
classes of multipartite entangled states. We first investi- 
gate the scaling properties of disentanglement of mixed-state 
generalizations of GHZ states, namely, generalized GHZ- 
diagonal states. For systems subject to natural decoherence 
models, analytical upper bounds can be derived for their en- 
tanglement and for the associated scaling behavior with the 
system size. Our findings are valid for any convex entangle- 
ment measure. In addition, we numerically study the entan- 
glement dynamics of randomly generated states, quantified 
in terms of the respective states' negativity. Random sam- 
ples of pure initial states do not abide to the above-mentioned 
bounds. Our numerical data suggest that the latter discrepancy 
increases with the number of system constituents, which gives 
evidence of the exponential fragility of GHZ entanglement to- 
wards decoherence not being a generic feature. 

The paper is organized as follows. In Sec. [TTJ we introduce 
our notation and define the generalizations of the GHZ states 
to be studied later. Sec. [Hi] defines our environment models. 
The analytical upper bounds for the entanglement, together 



with the resulting scaling behavior, are derived in Sec. IV 
where the tightness of these bounds is also assessed. Sec. |V| 
compares the scaling of the properties of GHZ-entanglement 
with those of random pure states. Finally, Sec. VI summarizes 
our conclusions. 
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II. GENERALIZED GHZ AND GHZ DIAGONAL STATES 

Originally, the GHZ-state [2] was defined for three qubits as 
a superposition of the many particle states of all of them being 
in their respective basis state |0) and all of them being in the 
orthogonal state This notion straightforwardly extends to 
more qubits and gives the GHZ-state to be 

I^ghz) = -^(|000... 0) + 1111...1)). (1) 

Without imposing more details on the nature of the specific 
quantum system the basis states |0) and |1) are abstract and 
without any physical meaning. Thus, formally, the state 
(|010) + 1 101)) /a/2 for three qubits, e.g., follows the spirit 
of the GHZ construction in precisely the same manner. When 
considering all such states it is helpful to interpret a string 
of zeros and ones as the binary representation of a number, 
the number of the second term always being bit-wise inverted 
with respect to the first. A many particle basis state of N 
qubits can thus be labeled with numbers from (all zeros) to 
2 N — 1 (all ones). Our three qubit example would therewith 
read(|2) + |5))/\/2. If we furthermore allow for different am- 
plitudes of the vectors and, for later convenience, treat a phase 
difference of it separately (thus distinguishing states of even 
and odd parity), we arrive at the generalized GHZ states of N 
qubits, 

|^±(a,/3)>=a|*>±/3|fc>, (2) 

with an N-bit number k, < k < 2 N ~ 1, k the bit- wise 
flipped number, and complex amplitudes a and (3 such that 
H 2 + |/3| 2 = 1, anda,/3 ^ 0. 

An incoherent mixture of several generalized GHZ-states 
|^ijr(a, j3)) with different k and parity, but the same ampli- 
tudes a and j3, is a generalized GHZ diagonal state 0361 : 

2 N -l 

P= E (A^ + («,/3))(^(a,/3)| + 

k=0 

+ X-\i:-(a,(3))(^(a,(3)\). (3) 

The coefficients denote the respective probabilities with 
which the state appears in the mixture. Naturally, they are 
positive and sum up to one. 

III. NOISE MODELS 

Given the above class of initial states, we proceed now 
to describe three paradigmatic models of incoherent dynam- 
ics present in typical experimental settings. We assume that 
qubits do not interact with each other, neither directly, nor 
indirectly through their baths, i.e., they feel independent and 
local environments. This is a reasonable approach for qubits 
being located in different laboratories, or sufficiently well iso- 
lated from each other. Furthermore, we assume that the qubit- 
bath interaction is identical for all of them. 



We describe the state evolution by means of a map (or chan- 
nel) such that a single qubit's initial state pi is mapped onto its 
final state by virtue of pt = £{pi). Since the mapping is be- 
tween quantum states, it needs to be completely positive and 
trace preserving. Every such map allows for an operator sum 
(or Kraus) representation 1 37 1 

£(p i )=Y^K j p i K}, (4) 

3 

with Kraus operators Kj fulfilling ^ . KjKj = 1. The initial 
joint state p of N qubits evolves thus according to the TV-fold 
tensor product of the individual maps: 

p = K{p)=£®£®---®£(p). (5) 
1 v / 

jv 

Specifically, the channels we consider ll37ll are: the depolar- 
izing channel (that occurs, e.g., in spin scattering), the dephas- 
ing (or phase damping) channel (typical of elastic collisional 
interactions), and a thermal bath at arbitrary temperature (pro- 
vided, e.g., by a thermal radiation background). 

The depolarizing channel (for which we will use the sub- 
script D in the sequel) describes the situation in which the en- 
vironment isotropically destroys the information on a qubit's 
state and thus steers it into a maximally mixed state 1/2. The 
characteristic quantity describing the dynamics is the proba- 
bility p of finding the state completely depolarized. The single 
qubit's initial state pi thus evolves to £u(pi) = p;(l — p) + 
tp/2. The corresponding Kraus representation allows for a 
convenient form in terms of the Pauli matrices: 

3 

£d(Pi) = E S J a 0Pi a 3 ! ( 6 ) 
3=0 

where s = 1 — 3p/4, si = s 2 = s 3 = p/A, a = 1, and 
a i, (72, and (T3 are the three familiar Pauli operators associated 
with the respective qubit. 

As a second type of noise, we discuss the phase damping 
channel (subscript PD). It represents a situation when quan- 
tum coherence is lost without any population or excitation ex- 
change. If the probability of complete phase loss is p, a Kraus 
representation reads 

£m(pi) = (1 -p)Pi+ 

+ p(|0)(0|p i |0)(0| + |l)(l|p i |l)(l|) (7) 

with Kraus operators 

M = iyr i P ) ^i = |0)(0|^ ) M 2 = |l)(l|Vp. (8) 

Finally, the third type of environment to be dealt with here 
is a thermal bath (subscript T). In this scenario, the qubit's 
basis states |0) and |1) function as ground and excited state, 
respectively, in order to exchange excitations with the bath. 
In the Born-Markov approximation ll38l [39l this yields the 
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Kraus operators 

^ = v^(|o><o| + |i)(i|Vw), 

K ^\fWi m ^ p > (9) 

K2 = sJ (|0)(0|VT^+ , and 

where the average excitation n of the bath modes induces an 
energy exchange from the bath to the system, and vice versa. 
For zero temperature, i.e., h — 0, this reduces to the ampli- 
tude damping channel, where system excitations irreversibly 
dissipate into the bath, triggered by the vacuum fluctuations 
of the bath mode, with a rate 7. The opposite case, at infinite 
temperature, is established by taking the limit n — > 00 and 
7^0, with 7J7 = r = const, and models a purely diffu- 
sive environment with diffusion constant T. In this case, p is 
the probability for an excitation exchange between system and 
bath. 

We stress that in all three cases, the models presented can 
encompass many different dynamics, depending on how one 
relates p to time t. For instance, the thermal bath at zero tem- 
perature can model an atom interacting with the free electro- 
magnetic field if we take p = 1 — cxp(— "ft/2), or with the 
field in a cavity by taking p = sin 2 (wi/2), in which Rabi 
oscillations (of vacuum Rabi frequency lo) can take place. 



IV. BOUNDS ON THE ROBUSTNESS OF 
ENTANGLEMENT 

With the detailed description of the states and the dynamics 
to be scrutinized at hand we can assess the robustness of the 
system's entanglement by bounding it from above as it evolves 
according to such open dynamics. We only assume two basic 
properties of the entanglement quantifier E and thus provide 
results for any such entanglement quantifier. We merely re- 
quire E to vanish for separable states, and not to increase 
when probabilistically mixing two states a and uj, i.e., it needs 
to be convex: 

E[pcr + (1 - p)uj] < pE(cr) + (1 - p)E(uj) (10) 

for any p S [0, 1]. 

The basic idea behind the derivation of the upper bounds 
is to decompose all studied states as convex combinations of 
an entangled and a separable part, i.e., p = p, en tp ea i + (1 — 
Ment)p.sep, an d to use tne properties listed just above to bound 
the entanglement evolution as E(p) < p en iE(p en i)- 

In what follows we identify such decompositions for initial 
states of the general GHZ or general GHZ diagonal type, un- 
der the influence of the paradigmatic environments considered 
in the previous section. 



A. Depolarization 

We begin with the case of open system dynamics that in- 
dependently depolarizes each qubit. For this type of decoher- 
ence process, we prove now the following: 

( i) The entanglement of an initially generalized GHZ diag- 
onal state p subject to local depolarizing Ad is bounded 
from above as 

E(A D (p))<(l~p) N E(p). 

To prove this statement, let us first begin by considering a 
single generalized GHZ state, namely the particular case po = 

The evolved density matrix A D ( / o ) is straightforwardly ob- 
tained by the iV-fold application of the depolarizing chan- 
nel |6]) of one qubit according to Eq. Q. In the resulting 
iV-fold operator sum we first focus on one of the sums for 
a single qubit. There the two terms for <j\ and 02 yield contri- 
butions of equal weight (s\ — s^) with different signs in their 
coherences, because o\ does a bit-flip whereas 02 results in 
a combined bit- and phase-flip (parity change). Altogether, 
this cancels the coherences and thus results in a diagonal and 
hence separable contribution. Thus, again considering all the 
sums, the only terms not immediately causing separable con- 
tributions are the applications of only the identity (cto) or the 
parity changing operator a%. 

Therefore, after the application of channel A D , we iden- 
tify three contributions to the final state. First, there is 
the unchanged state \^pQ (a, (3))(ipQ (a, (3)\ resulting from the 
application of only identities or an even number of parity 
changes. Second, there is the parity changed counterpart 
\ipQ(a, (3)}(ipQ(a, (3) | that is generated by an odd number of 
parity changes. Lastly, there is a separable contribution origi- 
nating from the application of at least one cr\ or ui to one of 
the qubits: 

A D (po) =A + |^ +(a,/3))(^ + (a,/3)| 
+A_|^o(«./Wo («./?) I 

+ A se pA i ep • (11) 

Here, the first two terms are of opposite parity. Their coher- 
ences partially cancel when being summed, yielding another 
diagonal (separable) contribution. Their difference is what de- 
termines the remaining entangled contribution in the decom- 
position. Its careful evaluation (see Appendix [A| yields 

AdOq) =(l-p) N Po + [l-(l- p )W]^ ep , (12 ) 

This result does not depend on the fact that it was derived 
using \ip^(a, (5)) with k — as initial state and holds for dif- 
ferent k as well. Also, since the application of channels is 
linear this result extends immediately to any convex combi- 
nation such as the generalized GHZ diagonal states given in 
Eq. (j3J. The convexity property ( fT0| ) then yields the desired 
result. 

As a corollary of ( i), we have: 
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( ii) For the special case of two qubits (N = 2), the bound 
in ( i) holds for all initial two-qubit states (not only the 
general GHZ diagonal ones). 

Consider first the pure-state case. Any pure two-qubit state 
|\E') can be expressed in local product-state bases such that it 
is |*) = a|00) + /3|11) = \ip£(a,0)), with real and posi- 
tive a and (3 (37 1. Then, for any pure two-qubit state, there is 
always a local basis (the Schmidt basis) in which it is a gen- 
eralized GHZ state. Now, since is basis-independent - it 
shrinks the Bloch sphere without distinction of any particular 
direction -, the local basis to which all four Kraus operators 
in |6]) make reference can be any, in particular that in which 
\ < S>) is a generalized GHZ state. By employing (i), we have 
that (ii) holds for all pure two-qubit states. 

For the mixed state case, we consider any two-qubit state 
p in its convex roof |40| decomposition p = XL Pnl^nK^nl 
(for which the entanglement and the average entanglement co- 
incide) and repeat the same reasonings used before. Linear- 
ity implies that A D (p) = E„PnA D (|*„)(*„|). Convexity 
implies that E[A D (p)} < E„Pn^[A D (|*„)(*„|)], with the 
latter being the average entanglement of Au(p). The latter is 
in turn smaller or equal than ^ n p n (l — fO^G^n))' for (i) 
holds for any pure two-qubit state, so in particular also for all 
of the |\& n ). Therefore, we have: 

E(A D (p)) < "£p n (l-p) 2 E(\* n )) 

= (l-p) 2 E(J2Pn\*n)) 

n 

= (l- P ) 2 E(p), (13) 

where from the first to the second line of ( fT3j ) we used the 
equivalence between the average entanglement and the entan- 
glement itself for the optimal decomposition. 

It is important to mention that, as was anticipated in the 
introduction, bound (ii) - mathematically expressed in ( ff3| l 
- is directly connected to a universal law discovered in 114 II . 
There, a universal bound on E[A(p)] was set for the particular 
case of E being the concurrence, but for £ any completely- 
positive map. Bound (ii) thus generalizes [41] in terms of the 
allowed entanglement quantifier E. It restricts, however, the 
environment coupling to the particular case of depolarization 
dynamics and yields a slightly weaker bound when evaluated 
for E being the concurrence. 

We also stress that bounds ( i) and ( ii) are optimal with re- 
spect to the class of states and entanglement quantifiers we 
deal with. This can be seen directly by showing an entan- 
glement quantifier saturating the bound for at least one state, 
which was done already in Ref. [31], where it was shown that 
the most resistant negativities of evolved GHZ-states tend to 
(1 — p) N times their initial value in the small p or large N 
limits. 



B. Dephasing 

For individual dephasing environments the results are sim- 
ilar to the previously obtained. In this case we prove the fol- 
lowing: 

( Hi) The entanglement of an initially generalized GHZ diag- 
onal state p subject to local dephasing Apo is bounded 
from above as 

E(A PD (p)) < (l-p) N E(p). (14) 

The strategy to show this is the same as in (i), but now 
with the dephasing channel |7]). In this fashion, we first 
identify the terms in A PD (p ) for a single general GHZ state 

Po = IV'o"( a i P))(' l Po( a ^ P)\ triat sum U P to some separa- 
ble state. We notice then that only the first Kraus operator 
of £pd (proportional to the identity operator) yields a non(- 
necessarily)-separable state upon application on (a, f3)). 
All other terms contain at least one diagonal projector, |0)(0| 
or |1)(1|, that eliminates the coherences of \ipq(cc, /?)) and 
takes it to some diagonal and hence separable matrix. There- 
fore, a decomposition of the final state is 

Apd(po) - (1 - p) N Po + [1 ~ (1 ~ p) N ]p sep , (15) 

where again p sep is some fully separable density operator; and 
where the procedure just used can be applied once more to 
any generalized GHZ-state \ipf(a, (3)). The rest of the proof 
follows thus as in ( i) from the linearity of the channel and the 
convexity property of the entanglement quantifier. 

We should note that, in contrast to the depolarizing channel, 
the dephasing channel is indeed basis-dependent - it does not 
commute with all local unitaries. Even though any pure two- 
qubit state is of generalized GHZ form in its Schmidt basis, 
this basis is not necessarily the one in which the dephasing 
channel is defined. Thus (Hi) does not generalize to all two- 
qubit states as was the case with (i). We also stress that again 
bound ( Hi) is tight, as follows from a similar argument as for 
bounds (i) and (ii). 



C. A thermal bath 

For the last example of environment models treated in this 
work we prove the following: 

(iv) The entanglement of any initially generalized GHZ- 
state pk = \ipk(ot,f))} (ipk (ot,[3)\ subject to local ther- 
mal baths At with an average excitation n in the bath 
modes is bounded from above as 
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E[Aj( Pk 



< 



N- 



1 - 



2n + l 



1 - 



1 



2n + 1 



p) +\P\ 2 1 



2n + 1 



1 - 



1 



AT-k 



2n + l 



where k is the number of ones in the binary string k, i.e. the 
number of excitations in the state \k), and E max is the maximal 
value that the entanglement quantify er E can take. 

As before, we start by considering the final state in the op- 
erator sum representation of the channel acting on the initial 
state and identify terms that cancel to separable contributions 
to the state. Similarly to (Hi) we find that only terms with 
the two Kraus operators Kq or out of all four of d9}, ap- 
plied to all of the qubits constitute the contributions that ren- 
der a non(-necessarily)-separable state. That is, Aj(p k ) — 
T, ju ... jN =o,2 K h ®...® K jN p k Kj t <g) . . . <8> K] N + A sep p sep . 



For brevity we refer with A en t/5ent to the first, in general non- 
separable term. By virtue of the entanglement quantifier's 
convexity the final entanglement is bounded from above, as 
in the previous cases, by the probability A en t of the state con- 
tribution pent in the first term in the convex sum, times its en- 
tanglement, and, trivially, by the maximal entanglement: 

E[A T (p k )} < X eat E(p ent ) < A ent £ max . (16) 

The probability A ent is given by the trace of the first term. It is 



Aent — Tr 



Tr 



K]K n ®---®K] N K 3N p 

ji,...jV=0,2 

n ( i i - .. " ,A io> <oi + ( 1 - • \p ) 11 >ci : i p, 



2n+l 



N 



= (<*,0) in 



2n+ 1 



" p)|0)(0|+fl-^rp]|i)(l| 



\a\ 2 I 1 



2n+ 1 



2n+ 1 

N- 



2n + I 



1 _ fl+ 1 
~ 2n + l P 



11 \ I n + 1 
P) 1 ~ n - , T P 



2n + l 



2n+ 1 



N-h 



(17) 



where the explicit definitions of K a and K 2 in Q were used 
in the second equality. 

The scaling factor in the bound of (iv) includes k as well as 
the amplitudes a and (3, and thus still incorporates details of 
the initial states. In order to arrive at a scaling factor that only 
includes details of the thermal environment and, of course, 
the number of qubits, we maximize the expression over all 
states \%jrj:(a, /?)) which we are focusing on. The solution is 
reached, for instance, in the limit \a\ 2 — > 1 and k — 0. This 
allows us to get a larger, but state independent, scaling factor 
such that we can rephrase (iv): 

(v) The entanglement of any initially generalized GHZ- 
state pk — \tp^(a, P))(tp^(a, P) \ subject to local ther- 
mal baths At with an average excitation n in the bath 
modes is bounded from above as 

/ \ N 

E[A T ( Pk )] < (l-^rP) ^max. (18) 

It is worth stressing that for this kind of noise model we 
are restricted to the (pure) generalized GHZ-states. Also, no- 



tice that for different temperatures, and thus different mean 
numbers of bath excitations n, bounds (iv) and (v) range from 
E[A T (p k )] < [|cf (l-p) K + |/3| 2 (l-p) w - K ]£ max forbound 
(iv) and a trivial bound for (v), E\A^(p k y\ < E mRX , in the 
purely-dissipative, zero-temperature limit with h = (ampli- 
tude damping channel), to i?[A T (p fe )] < (1 — ^) N E max for 
both bounds in the purely-diffusive, infinite-temperature limit 

71 — > CO. 



V. BEYOND THE GENERALIZED GHZ STATES: 
COMPARISON WITH RANDOM STATES 

The simplicity of the obtained bounds raises the question of 
how dependent they are on the specific choice of initial states. 
Some of the above bounds are in this sense more general than 
others, but all feature an exponential scaling in the number N 
of qubits, which renders their entanglement more fragile for 
increasingly many qubits. 

In order to verify up to what extent the bounds are represen- 
tative to a larger class of states than only to the ones they were 
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Figure 1: (Color online) Normalized negativities of most balanced bipartitions as a function of p for systems of different size N undergoing 
individual depolarization. The red solid line corresponds to the bound (1— p) , the dotted line to the balanced GHZ-state \if)g (l/y/2, 1/V2)), 
which always lies below the bound (1 — p) and approaches it as N grows, and the blue crosses to the average over 10000 initially random pure 
states distributed uniformly over the system Hilbert space. The gray shadings are the distributions of the normalized negativities around their 
mean values (bin population reflected by saturation, i.e. dark is the maximum and white the minimum). For TV = 6 the average normalized 
negativity no longer lies below the bound (1 — p) . See text. 
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Figure 2: (Color Online) Same as Fig. [I] but for the least-balanced bipartitions. In this case the violation of the (1 — p) N scaling law already 
appears at N = 5. See text. 



actually derived for, we consider initial random pure states 
uniformly distributed over the entire Hilbert space j42j . As a 
simple, convex entanglement quantifier we choose negativity 
evaluated for a bipartition of the N subsystems. For systems 
with varying number of constituents, we generate a sample of 
10000 initial states, evolve it under the three considered inco- 
herent dynamics, and calculate the negativity as a function of 
the probability p of an incoherent event. 

In general, we observe that such states violate the above 
bounds, more drastically as we increase the number of parties 
N. Figs. [T] and [2] show examples of this violation. There, the 
average negativities Neg(p) (normalized to their initial value 
Neg(O)) of a sample of 10000 initial pure states undergoing 
individual depolarization is plotted for systems of different 
size N (blue crosses). For comparison, they also contain the 



bound (1 — p) N (red solid line) and the specific case of the 
balanced GHZ state \ip£ (l/y/2, I/a/2)) (dotted line). Nec- 
essarily, the latter always lies below the bound (1 — p) N and 
approaches it as N grows. The average normalized negativity, 
however, violates the (1 —p) N exponential-decay scaling law, 
in particular for larger N. The gray shadings along the vertical 
direction represent the histogram of the samples' normalized 
negativities in gray-scale (bin population reflected by satura- 
tion, i.e. black the maximum and white the minimum) ||431 . 
Fig. [T] presents the data for negativity evaluated for the most 
balanced partition of N/2 versus N/2, where the violation is 
apparent for 6 qubits. Contrary, Fig. [2] shows the least bal- 
anced partition of 1 versus N — 1 qubits, where a violation 
appears for 5 qubits. 

With the observation that initial random states violate the 
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N 


Sample size 


2 to 7 


10000 


8, 9, 10 


5000 


11 


1560 


12 


100 


13 


10 


14 


1 



Table I: Sample size used in the numerical calculation of Fig. [3] 

bounds, we can consider an extension of these bounds in such 
a fashion that these states are included as well without losing 
the scaling behaviour, i.e., an entanglement upper bound that 
decreases with increasing N for all the states. Fig.[3]suggests 
that this is not possible - the mean value of the normalized 
negativity, of the least balanced bipartition, at a given evolu- 
tion step (dephasing at p = 0.3) grows with increasing N. In 
this plot, the full line represents the the bound (1 — p) N , while 
the dots stand for the numerical values obtained by a similar 
sampling as explained above. However, due to the big com- 
putational effort of sampling many states of large N, we rely 
on the negativity concentration effect for large dimensional 
systems under incoherent dynamics ll43l l44l |451 in order to 
pick just few typical states. The sample sizes are displayed in 
Table |U 

A persistently increasing mean negativity, which we con- 
firmed numerically up to N = 14 qubits, can thus not be 
bounded from above by any nontrivially decreasing function. 
This also implies that the exponential fragility of GHZ-type 
states is not typical, as it is not a decisive feature present in 
the ensemble of randomly chosen initial states. 

The increasing entanglement robustness for random states 
is observed for all channels here scrutinized. However, it only 
takes place for highly unbalanced partitions, where the en- 
tanglement between few qubits with the remaining system is 
probed. For more balanced partitions, although the bounds 
presented in the previous section are violated (as shown in 
Fig. |TJ, we observe a stronger entanglement decay of random 
states for increasing number of subsystems. Thus, if an im- 
proved bound exists which applies also to random states, it 
must take into account the relative partition sizes in a many- 
body state 



Neg(0.3) 
Neg(0) 0.3 




2 3 4 5 6 7 8 9 10 11 12 13 14 
N 



Figure 3: (Color Online) Mean normalized negativity for the most 
unbalanced partition as a function of the system size, for a fixed value 
of the probability p — 0.3, and random initial states undergoing the 
action of the dephasing channel (circles, joined by straight lines as a 
guide to the eye). Bound (Hi) is plotted as a solid line. 



present conclusions are not only valid for bipartite but also for 
any genuine multipartite entanglement the system can present. 
Moreover, most of the entanglement quantifiers aiming at the 
usefulness of a given quantum state for certain applications 
that involve local operations fall into the category of quanti- 
fiers adopted in the present work |46l . In this way, our results 
bound the GHZ-states' usefulness for most of quantum com- 
municational and computational tasks. 
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VI. CONCLUSIONS 

In this paper we have considered generalized GHZ and 
GHZ-diagonal states evolving under the action of natural de- 
coherence processes. For these states we have derived analyt- 
ical upper bounds to the entanglement decay and shown that 
the fragility of their entanglement increases exponentially as 
the size of the system grows and independent of relative par- 
tition sizes. A comparison with random pure initial states un- 
der the same dynamics manifests the exponential fragility of 
entanglement towards decoherence to be a distinct feature of 
GHZ-type states. 

We stress once more that the class of entanglement quanti- 
fiers considered in the calculations is very general. Thus, the 



Appendix A: ENTANGLED CONTRIBUTION IN THE 
OPERATOR SUM REPRESENTATION OF DEPOLARIZING 

In the proof of ( i), the application of the depolarizing chan- 
nel Ad to the initial state po = \i>o ( a i @)) (^o ( a > @)\ yields, 
in the operator sum representation using |6]), two contribu- 
tions. The application of operators o\ or a 2 to any of the 
qubits gives a separable contribution, whereas a non-separable 
contribution can exclusively appear for the application of only 
<7o = 1 or 03 to all the qubits: 

s J'i • ■ • s m a ji®- ■ -®^ 3 NP^ n ®- ■ -®o~ jN (Al) 

Jii — j'jv=0,3 
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with respective prefactors sq — 1 — 3p/4 and S3 = p/4. An 
even number of parity changes due to 03 (and identity oper- 
ators everywhere else) leaves the state effectively unchanged, 
|^|q (a, (3))(i^q (a, /3)\, whereas an odd number changes the 
parity, \i/jq (a, 0)){i/jq (a, (3)\. The corresponding prefactors 
are then accordingly given by the sum over all even and odd 
powers of s 3 and s > respectively, weighted with their oc- 
curence: 

M=0,even x 7 v 7 
A/=l,odd v 7 v 7 



Their difference is what determines the contribution to the re- 
maining coherences, by virtue of the binomial theorem: 




(A4) 



= (1-P) . (A5) 
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